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Abstract 

The role of orbital ordering on metal-insulator transition in transition-metal oxides is inves- 
tigated by the cluster self-consistent field approach in the strong correlation regime. A clear 
dependence of the insulating gap of single-particle excitation spectra on the orbital order parameter 
is found. The thermal fluctuation drives the orbital order-disorder transition, diminishes the gap 
and leads to the metal-insulator transition. The unusual temperature dependence of the orbital 
polarization in the orbital insulator is also manifested in the resonant x-ray scattering intensity. 
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I. INTRODUCTION 



Metal-insulator transition (MIT) in strongly correlated transition-metal oxides (TMO) is 
one of the central problems in condensed matter physics, it has attracted extensive attention 
in recent decades since it has been found that high temperature superconductivity, colossal 
magnetoresistance and many other phenomena occur in the vicinity of MIT 1 ' 2 ). These 
strongly correlated electronic systems exhibit very complicated and rich phase diagrams with 
temperature, doping, pressure and magnetic field 3 ' 4 . In these compounds, the temperature 
induced MIT in V 2 3 and manganites is especially interesting, because the MIT temperature 
T M is much smaller than the insulating gap (A) in transport, for example in V 2 3 , T M =154 
K<<A = 0.6 eV 5 \ Obviously, the thermal fluctuation is not the driven force of the MIT. 
It has recently realized that the orbital degree of freedom and orbital ordering (00) play 
important roles in the groundstate (GS) properties of these TMO. 

The 00 was first proposed to explain the complicated magnetic structures by Kugel et al. 
6 ) for KCuF 3 and by Castellani et al. 7 for V2O3. Due to the strong correlation between 3d 
electrons and the large anisotropy of the 3d wavef unctions, the orbital degree of freedom and 
the 00 affect many electronic and magnetic properties of these strongly correlated systems 
4 ' 6 ~ 9 \ and have been extensively studied in colossal magnetoresistive manganites, vanadium 
oxides and many other TMO [for example, see Ref.10]. Experimentally orbital order- disorder 
transition is usually accompanied by MIT, such as in prototype MIT compound V2O3, an 
obvious insulator to metal transition occurs at Tm=Tc n \ here Tc is the Curie- Weiss tem- 
perature. It is believed that in V 2 3 , the 00 transition (00T), the MIT and the magnetic 
transition occur simultaneously, i.e. T M = T OQ = T c 12 \ We also notice that in lightly 
doped Lao.88<SVo.i2-^i03 the ferromagnetic insulator to ferromagnetic metal transiting at 
Tm is identified as an 00T by Endoh et al 13 \ or, Tm = Too, ruling out the significant 
role of cooperative Jahn- Teller (JT) distortion on the GS. The orbital phase transition near 
the vicinity of MIT is also reported in Lai- x Ca x Mn0 3 (x « 0.2) 14 ). These experiments 
clearly establish a close relationship between 00T and MIT. Theoretically Castellani et al. 
first suggested the correlation between 00T and MIT in V 2 3 7 \ they proposed that the 
MIT associated with the 00 was driven by the variation of the entropy in the presence of 
long-range magnetic and orbital orders. Khaliullin et al. 15 attributed this correlation to 
the formation of orbital polaron. Nevertheless, it is not well understood theoretically how 
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the 00 insulator evolves to MIT with lifting temperature. 

In this paper, we demonstrate that the strong orbital correlation in TMO with orbital 
degenerate 3d electrons leads the GS to be an orbital insulator, and the orbital order- disorder 
driven MIT. Starting from a twofold-degenerate spin-orbital interacting model, and using 
the cluster self-consistent field (cluster-SCF) approach developed recently, we first determine 
the 00 GS; and then show that the insulating gap of the single-particle spectrum opens as 
the long-range 00 establishes. With increasing temperature, the thermal fluctuation drives 
OOT, also the energy gap of the single-particle excitation vanishes at T 00 , indicating the 
transition from insulator to metal. The resonant X-ray scattering (RXS) intensity diminishes 
to zero near the critical temperature Too- The rest of this paper is organized as follows: in 
Sec. II we describe the effective model Hamiltonian of an orbital insulator and the cluster- 
SCF approach; then we present the magnetic and orbital structures in GS, the T-dependent 
00 parameters and the gap of single-particle excitation spectra in Sec. Ill; the variation 
of the temperature and the azimuthal angle dependence of the RSX intensity are given in 
Sec. IV, and the last section is devoted to the remarks and summary. 

II. MODEL HAMILTONIAN AND METHOD 

In many perovskite transition-metal oxides under the octahedral crystalline field, the five- 
fold degenerate 3d orbits of the transition metal ions split into lower t 2g and higher E g orbits. 
For clarification and simplification we consider such an ideal cubic TMO system that the t 2g 
orbits are filled and contribute no spin, and the twofold degenerate E g orbits are occupied by 
one hole or one electron, corresponding to the electron configuration of 3d 9 in KCuF 3 or 3d 7 
in LaNi03. Such a system is spin-1/2 and twofold orbital degenerate, or the orbital pseu- 
dospin t—1/2. We denote the two E g orbits as |l)=| e gi)H^32 2 -r 2 ) an d \^)—\ e g2)—\d x 2 -y 2 )- 
The major low-energy physics of the system is described by the twofold-degenerate Hubbard 
model with strong Coulomb interaction 6 ' 7 \ In the second-order perturbation approximation 
the E g electrons interact with each other through the low-energy superexchange coupling, 
which is expressed as: 




l=x,y,z 




(2.1) 
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where the operator l\ = cos (2nmi/3) rf — sin (27rm;/3) rf, the index 1, I = x, y 
or z, denotes the direction of a bond; connects site i and its nearest-neighbor 

site j along the I direction, and (m x ,m y ,m z )— (1,2,3). t z and t x are the Pauli ma- 
trix, r z = | represents the orbital polarization in the state |1) and r z = — | 
the orbital polarization in the state |2). Thus the polarization degree of the or- 
bit, (r), is called the orbitalization. The constants J\, J2, J3 and J4 are the su- 
perexchange interactions, and are read: Ji = 8t 2 [U/(U 2 — Jjj) — J H /(U 2 — J%)], J2 = 
16t 2 [1/(17! + J H ) + 1/(U + J H )\, J 3 = 32t 2 [UJiUl - J 2 H ) - J H /(U 2 - J 2 H )}, and J 4 = 
8t 2 [(Ui + 2J H )/{Ul - J 2 H ) + J H /(U 2 - J 2 H )], with U = Ux + 2J H , here At is the hopping 
integral between the |2) orbits along the z direction. U and U\ are the intra- and inter- 
orbital Coulomb interactions, and Jh is the Hund's rule coupling. In this paper we adopt 
t= 0.1 eV and J H = 0.9 eV. 

Clearly, the fourth term in Eq.(l) is an orbital frustration interaction: while the exchange 
along the z direction stabilizes the alternating \3z 2 — r 2 ) and \x 2 — y 2 ) configuration, the 
equivalent coupling along the x or the y directions favors other orbital pair configuration, 
i.e. 1 3a; 2 — r 2 ) and \y 2 — z 2 ) or \?>y 2 — r 2 ) and \z 2 — x 2 ). The second term in Eq.(l) is 
a "magnetic field" for the orbital pseudospin, i.e. the orbital field. Due to the relation 
If + If + If = 0, the orbital field favors a peculiar orbital polarization and suppresses the or- 
bital quantum frustration. However, if the spin correlations (sjSj) are identical along the x, y 
and ^-directions, such as in the ferromagnetic (FM) or the Neel antiferromagnetic (G-type 
AFM) ordered structure, the orbital field vanishes. Then the third spin-orbital interaction, 
which are equal along the x, the y and the ^directions, together with the fourth term leads to 
the strong quantum fluctuations and results in an orbital disordered ground state 9 . Fortu- 
nately, the conventional JT instabiity lowers the cubic symmetry to the tetragonal symmetry 
through the JT distortion and lifts the orbital degeneracy. The tetragonal crystalline field, 
H z = E z J2i T ii is & n additional orbital field to break the symmetry of orbital space, and 
favors a peculiar orbital configuration. Thus the full Hamiltonian including the crystalline 
field term is read: 

H = H SE + H Z . (2.2) 

We will discuss the cooperative JT effect on the 00 in the Sect.V. 

It is still a difficult task to treat the spin and orbital correlations and fluctuations, and to 
find the GS of Eq.(2) with high accuracy. To study the GS properties of these TMO systems, 
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we recently developed the cluster self-consistent field (cluster-SCF) approach to study the 
00 properties and the evolution of GS properties with interaction parameters in V2O3 16 . 
This approach combines the exact diagonalization and the self-consistent field techniques, 
the main idea is briefly addressed as follows: firstly, choose a proper cluster, usually the unit 
cell of the TMO compound. Secondly, divide the interactions in Eq.(2) into two types: the 
internal bond interactions, H,j, between the sites Rj and Rj in the cluster, and the external 
bond interactions, H^, between the environment site R/ and the site Rj in the cluster. The 
periodic condition enforces the electronic states at the environmental site Ri to be identical 
to the corresponding atom inside the cluster. After diagonalizing the internal interaction 
Hjj of the cluster, the interactions of the environment site R; in acting on the cluster 
are obtained as a initial value of the self-consistent field (SCF): h—Tr^pnH'^), here pu is 
the reduced density matrix of the bond (il) after tracing over all of the other sites in the 
cluster. Thirdly, diagonalize the cluster Hamiltonian H^- in the presence of the SCF h[. 
Iteration is then performed until the orbital correlation functions and the SCF h- converge. 
The GS properties of the system and the effective field of the surrounding atoms applied 
on the inside atoms of the cluster are thus obtained. One of the advantages of the present 
approach is that the short-range correlation of the spins and the orbits, which is neglected 
in the conventional self-consistent mean- field method, is taken into account. 

To check the validity and convergence of the cluster-SCF approach, we apply it on the 
simple spin- 1/2 AFM Heisenberg model. The Hamiltonian is read: 



where Si is the operator for the i — th spin and J > 0. In our numerical calculation code, the 
convergence of the GS energy is considered as the truncation condition of iterations, when 
the relative error of the GS energy between two iterations is within 10~ 8 , the GS energy is 
converged. The net nearest-neighbor spin correlation function S(si, Sj)—{si ■ Sj) — (s*j) ■ (sj) 
and the spin polarization (s z ) are used cooperatively to determine the magnetic structure. 
We find S(si, Sj)—— 0.1115 along the three axes, and the spin polarizations of lattice clearly 
falls into two distinct kinds: (s z )=— 0.4683, and (s 2 )=0.4683. Obviously the GS is the Neel 
AFM structure, and the averaged spin deviation of the z component from 1/2 at each site 
is 0.0317, smaller than Anderson's zero-point fluctuation result 0.078 17 \ This difference is 
attributed to the lack of the long-range quantum fluctuation of the spins in our approach, 




(2.3) 



5 



the present cluster-SCF approach includes only the short-range quantum fluctuations. This 
result clearly confirmed the efficiency and validity of our approach. 

III. ORBITAL TRANSITION INDUCED METAL-INSULATOR TRANSITION 

In this section we first utilize the cluster-SCF method to find the magnetic and orbital GS 
of Eq.(2), then explore the evolutions of the ordered orbital parameter and the energy gap 
of single-particle excitation spectra with increasing temperature in the strong correlation 
limit. 

A. Magnetic and orbital structures in ground state 

In many three-dimensional TMO the magnetic structure s are not difficult to determine 
by the neutron scattering and other experimental techniques, in the following in accordance 
with the spin configuration in the experiments, we treat the large local spins as semiclassical 
to save computation time, and focus on the orbital fluctuation and the orbital GS for various 
magnetic orders. To determine the most stable magnetic and orbital GS structure of Eq.(2), 
we choose a cubic 8-site cluster. In the cluster, each site has three internal bonds and three 
external bonds. The total energy of the cluster as a function of the ratio U / Jh for various 
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FIG. 1: Dependence of the total energy on the ratio U/Jh for different magnetic structures. CO 
and GO denote the C-type and Neel antiferro-orbital ordered structures, respectively. Theoretical 
parameters: Jh = 0.9eV and t = O.leC, E z = —bmeV. 
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FIG. 2: Dependence of GS energy on the ratio U/Jh in the classical approximation. A-AFO 
denotes the A-type antiferro-orbital ordered structure. Theoretical parameters are the same to 
Fig.l 

magnetic structures is shown in Fig.l. It is found that for small U/Jh, the GS is FM, which 
agrees with the mean-field results of Roth 18 and Cyrot and Lyon-Caen's 19 ; while the GS is 
A-type AFM (A-AFM) for large U / Jh, he. AFM coupling along the c-axis and FM couplings 
in the a6-plane. In contrast, the classical approximation, which treats both the spin and the 
orbital operators as classical, shows that the GS is A-type AFM structure (A-AFM) over 
very wide U/Jh range, as seen in Fig. 2. This classical result is in agreement with Kugel and 
Khomskii's result 6 >. As we will show later, after taking into account the quantum effect, 
the orbital occupation and polarization are slightly different from the classical results. 

To find the different stable magnetic structure sand the corresponding orbital con- 
figurations in different U / Jh range, we take the net orbital-orbital correlation function, 
Cij=(ri ■ Tj) — (r/) ■ (tj), together with the orbitalization (r) to determine the orbital GS. 
In the FM phase, the Ising-like orbital-orbital interactions, which are identical for 

I = x,y and z, as we analyzed in Sec. II, favors the orbital liquid GS. The tetragonal crys- 
talline field suppresses the quantum fluctuations, drives the electrons into the \3z 2 — r 2 ) orbit 
in E z < 0, and stabilizes the GS as G-type antiferro-orbital (AFO) ordering, which is most 
favorable of the the anisotropic distribution of the electron clouds. This agrees with the 
empirical Goodenough-Kanamori rules 20 \ Our numerical results further show that all of 
the orbital correlations, 5y along the x, y and z-directions, are AFO, thus each site is AFO 
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polarized with respect to its nearest-neighbor sites. One finds that the orbital occupations 
in the two sublattices are the combinations of |1) and |2): (v^l) + >/l9|2» /ViO, and 
(v^l|l) - >/l9|2)) /y/40, for U/J H =4.0 and E z = -5meV. 

With the increasing of U/Jh, the electronic superexchange interaction becomes small, in 
comparison with the crystalline field splitting, and more and more electrons occupy the \x 2 — 
y 2 ) orbit. The FM phase becomes unstable and transits to the AFM GS at (£7/J#) c =4.1. At 
U I J#=6.0, the orbital occupations in the two sublattices are approximately (v3|l) + |2)) /2 
and (v3|l) — |2)) /2, as shown in Fig. 3. The single-bond coupling energy, which is calculated 
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FIG. 3: Single-bond spin coupling (SC) strength (3a) and orbital coupling (OC) energy (3b) as 
the function of U / Jh for FM and A- AFM structures. F x r y g \ and &- x (y,z) denote the x (y, z) bonds 
in the FM and the A-AFM structures, respectively. 



for specific spin configuration via E-j=Tr s (pijHij), consists of the orbital-dependent spin 
coupling energy and the pure orbital coupling energy. It is found in Fig. 3a and Fig.3b that 
for small U / Jh, the FM structure is the most stable, and the quantum fluctuation of spins is 
small; On the contrary, the A-AFM structure is the most stable for large U/Jh, addressing 
the GS orbital configuration in typical OO compound KCuF^. In the following we focus 
the small U/J# case. More detail results for large U/Jh and the effect of crystalline field 
splitting on 00 compound KCuF^ will be presented in a further paper. 

The orbital phase transitions with the variation of the Coulomb interaction are typical 
quantum phase transitions. We find that the concurrence of two nearest-neighbor pseu- 
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dospins critically changes at the quantum transition point, detail result about how the 
entanglement of the orbital states evolves with the quantum phase transition will be pub- 
lished in the future. On the other hand, if we freeze the orbital configuration as the G-AFO 
ordered phase and to search the most stable spin configuration, we find that the magnetic 
structure transition from FM to A-AFM also occurs at the same critical value (U / J#) c =4.1, 
confirming the validity of the orbital GS. 

B. Temperature dependence of orbital order parameter 

After determining the orbital GS, we explore the influence of the thermal fluctuations 
on the spin and orbital order parameters. At finite temperature T, thermal fluctuations 
exciting spin waves and orbital waves weaken both the orbital and the spin orderings, and 
destroy the magnetic order at the Curie temperature Tc and the 00 at the orbital critical 
temperature Too- Since the spin order interplays with the 00, the reduction of sponta- 
neous magnetization with increasing temperature also softens the orbital interaction, and 
vice versa. In the following we present the temperature dependence of the magnetic order 
and the 00 parameters in small U/J H range (U/J H < 4.1) in the conventional mean-field 
approximation. In the FM and G-AFO phase, the mean-field Hamiltonian is approximated 
as: 

HmF = (^A(B) T iA(B) + Jl(B) T iA(B) + J sS'i) (3.1) 

i 

where the coefficients J^my Ja{by and J s depend on the magnetic order and the 00 
parameters, for example: 

J s = 3J X <**) + 3/2J 3 (s z ) «rl> (r*> + (rl) (7%)) (3.2) 

etc. In the FM and G-AFO phase, the thermal averages of the spontaneous magnetiza- 
tion and the orbital polarization, i.e. the orbitalization, satisfy the following self-consistent 
equations: 

«) = Tr{s z a exp(-pj saS :)}/Z s a (3.3) 
(r x A [z) ) = Tr{rf ) e* V (-l3jf\f ) )}/zf ) (3.4) 

where (3=l/k B T. Also, the self-consistent equations for the A-AFM and G-AFO phase at 
large U / Jh can be obtained, but are more complicated. The temperature dependence of the 
spin and 00 parameters for the system with small U / Jh is shown in Fig. 4. 
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At low temperature, the magnetization and the two components of the orbitalization are 
nearly saturated. At lifting temperature T>Tc/2, the magnetization decreases considerably, 
while the orbitalization (r) almost does not change. With the temperature approaching Tc, 
the magnetization rapidly falls to zero, and the magnetic transition is obviously the first 
order, which originates from the strong anisotropy of the spin-orbit coupling, in agreement 
with T. M. Rice's argument based the Landau phenomenal theory for multiple parameter 
orders 22 ^ and our renormalization group analysis to the critical indexes of the spin-orbital 
interaction models 23 \ Meanwhile the 00 parameter steeply diminishes a fraction at Tc, 
as observed in Fig. 4. With further increasing in temperature, the 00 parameter gradu- 
ally decreases and finally vanishes at T o, with the characters of the second-order phase 
transition. 

Due to the existence of orbital-orbital interaction in Hamiltonian (1), the long-range 
orbital correlation still exists after the magnetic order disappears. Therefore the orbital 
order-disorder transition temperature Too is higher than the magnetic transition tempera- 
ture Tc, as observed in Fig. 4. It is worthwhile noticing that the orbitalization components 
(t z ) and (t x ) display considerable discontinuity at the Curie temperature Tc, which indi- 
cates the strong correlation between the spin and the orbital degrees of freedom. The falling 
fraction of the orbitalization near Tc is about 12% with respect to the saturated magnitude. 
The ab initio electronic structure calculations based on the local density approximation with 



0.5 

2 0.4 
CD 

CD 

cfl 0.3 
i_ 
cfl 

a 

£ 0.2 
"O 

i— 
O 

0.1 



100 200 300 

T(K) 

FIG. 4: Temperature dependence of the magnetization s z , the orbital polarizations t z and t x . 
Theoretical parameters: U = 3.6eF, Jh = 0.9eF, t = O.leV '. 

correlation correction (LDA+U) have shown that in KCuF 3 the variation of the orbitaliza- 
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tion from the low-temperature AFM phase to the high-temperature paramagnetic phase is 
small, about 5% to 10% 24 \ Experiments, however, showed the variation of the orbitaliza- 
tion is about 50% or even 100% in KCuF 3 22 and in V 2 3 25 . This discrepancy may arise 
from two reasons: the first one is that the paramagnetic phase calculated by the LDA+U 
approach at zero temperature neglects the change of the entropy due to the variation of the 
spin and orbital configurations; the second one is that in the LDA+U approach and the 
present work the long-range orbital fluctuations are not taken into account. 

C. Energy gap of single-particle spectra 

In this subsection we study how the gap of single-particle energy spectra of the E g elec- 
trons depends on the evolution of the 00 parameter. In the OO phase, the long-range 
orbital order breaks the symmetry of the orbital space, and opens an insulating gap A in 
single-particle excitation spectra of the TMO, forming an orbital insulator. Such an insu- 
lator differs from the single-band Mott-Hubbard insulators in many aspects. One of the 
outstanding characters is that the insulating gap strongly depends on the OO parameter, 
hence closely on the temperature. In the TMO with quarter filling, the tight-bonding spec- 
trum of an Eg electron reads: 

1 dj2 + d\ 2 dj\^ + 36^2^^ 
+ AtJ2( d lidji + h.c.) (3.5) 

where dj creates a 3d electron at site % with orbital state 7; and ± denote the signs of the 
hopping integrals between two nearest-neighbor sites along the x and y direction, respec- 
tively. We ignore the spin index in the present FM GS. 

Strong on-site Coulomb interaction U is the main character in these ordered TMO, that 
U >> t prohibits the E g electrons from double occupation at the same site, thus the single 
occupation constraint is applied on the E g orbits. In the limit of large Coulomb interaction, 
the constraint of no double occupancy at site Ri is enforced by introducing auxiliary fermions 
26 \ fi 7 , an d bosons, bf, here f% creates a slaved fermion (electron) at site Ri with orbital 
state 7, and b\ creates a boson (hole) at site thus d\ ia —f} lc bi, and the single occupation 
condition is: J2a-y fijafho- + 1- m ^ ne FM and G-AFO ordered phase, the effective 
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low-energy Hamiltonian becomes: 

Heff = E[4 1 / fe t 1 /.i + 4 2 / fe U2 + el 2 (flf k2 

k 

+ fLfkl) + Jl(fk+Q,lfk2 + fk+Q,2fkl)} 

+ E A ^i^i + / fc V fc2 + & 2 -l) (3.6) 

k 

here Q=(vr, 7r, 7r), corresponding to the G-AFO order. The average of the boson occu- 
pation is approximated as a c-number: < b\bi >= x. The dispersion functions are: 
e^ 1 = 4rt(cos k x + cos k y +4cosk z ) + 3 (J 3 (s z ) 2 + J 4 ) ((r 2 ) + E z /2), ef = 4 v / 3a:t(- cos k x 
+ cosk y ), and e 22 = 12xt(cos k x + cosk y ) +3(J 3 (s z ) 2 + J 4 ) ((r 2 ) + E z /2). And the parameter 
Ji=— 3 (J3 (s z ) 2 + J4) (r x ). In the saddle-point approximation, the constraint constant A 
and the boson occupation are: 

<»,J>77' 

^ = 1 -^E((/i/«) + (/i/»)) ^ 

i 

respectively. Physically A shifts the energy level of the Fermi quasiparticles, and x gives 
rise to the hole concentration upon doping. In the compounds with quarter-filling, x — > 0. 
The single-quasiparticle energy spectra E k have four branches and strongly depend on the 
00 parameter. For the G-AFO ordered TMO the lower two subbands are filled, since each 
magnetic and orbital unit cell contains two electrons at quarter-filling, the minimum of the 
empty subband separates from the maximum of the filled subbands by a gap A. This result 
considerably differs from Kilian and Khaliullian's gapless charge excitation 15 ); in fact, in 
their study the modulation of the long-range 00 on the motion of the electrons was not 
taken into account. 

The thermal fluctuation weakens the 00 parameter, hence softens the single-particle 
excitation spectra. The T-dependence of the energy gap of the single-particle spectra is 
shown in Fig.5. At low T, the strong Coulomb interaction of the 3d electrons, the large 
spatial anisotropy of 3d orbits and the crystalline field splitting cooperatively lift the orbital 
degeneracy, and favor the ordered orbital phase. The broken of the orbital space symmetry 
leads to that an electron must cost at least A in energy to hop to the nearest-neighbor 
site. Therefore the threshold A gives rise to the insulating gap of the TMO as an orbital 
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FIG. 5: Temperature dependence of the energy gap of the single-particle spectra in the ferromag- 
netic and G-type antiferro-orbital ordered phase. A = — O.leV, and other parameters are the same 
to Fig.4. 

insulator. This energy gap manifests in the optical excitation and transport, and differs from 
the energy gap of the orbital wave excitation. Increasing temperature gradually destroys the 
long-range OO, and leads to the decrease of the gap of the separated subbands. With the 
further increasing of temperature to T = Too, the thermal fluctuations become so strong 
that the long-range OO and the OO parameter vanish. In this situation, the two lower 
subbands merge together with the two higher subbands, the insulating gap disappears, see 
Fig. 5, indicating the occurrence of MIT in TMO. Thus the energy gap of the orbital insulator 
crucially depends on the OO parameter, hence the temperature. 

The temperature dependence of the insulating gap A(T) turns out to be similar to that 
of the OO parameter, steeply diminishes a fraction at the Curie temperature Tc and com- 
pletely disappears at Too, see Fig.4 and Fig. 5. Therefore, the temperature driven MIT in 
TMO is essentially induced by the orbital order- disorder transition. Our numerical results 
demonstrate the significant discrepancy between the optical and transport gap A(T) and 
the MIT critical temperature T M . For the present system in Fig. 5, the gap of the energy 
spectra, usually reflected in optical and transport experiments, is A = 0.5 eV; while the 
MIT critical temperature Tm = Too ~ 270K = 0.023, which accounts for the considerable 
discrepancy between A and Tm in many orbital compounds. 
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IV. TEMPERATURE-DEPENDENT RESONANT X-RAY SCATTERING IN- 
TENSITIES 

Strong interplay between spin and orbital degrees of freedom also manifests In the op- 
tical excitation of orbital insulating TMO. Recently it was proposed to utilizing the RXS 
technique to measure the 00 phase by Murakami et al. 28 for manganites and Fabrizio 
et al. 29 for V2O3. The Is — 3c? K-edge RXS peak intensity is a useful signal to probe 
the 00 phase and the orbital order- disorder phase transitions n ' 27 . The polarization and 
azimuthal dependences of the RXS intensities provide the hidden information of the under- 
lying 00 parameters 27_30 ) ; and unveil the interplay between spins and orbits. Although 
the signal enhancement of the quadrupole Is — 3d scattering (E 2 ) is less than that of the 
electric dipole Is — 4p scattering (Ei), the observed E2 spectral line shape as a function of 
energy is more easy to be identified than the quite complicated line shapes associated with 
Ei process n > 28 > 29 ). In what follows we present the E 2 RXS intensity and its evolution with 
temperature based on the 00 phase obtained in the preceding section. 

In the FM and G-AFO phase, the ordered 3d orbits in the sublattices consist of two 
different orbital basis: |?/>i) = oki|1) + a 2 \2) and \ip 2 ) — — ck 2 j 2) , here the coefficients 

ai,a 2 are the functions of the interaction parameters. There exist two kinds of reflections: 
the fundamental reflection at (hkl) = (n x n y n z ) and the orbital superlattice reflection at (hkl) 
— { n x + \i n y + \i n z + §)• At the lattice symmetry-forbidden direction in which (hkl) are 
all odd, the orbital structure factor is written as: 



where the function / (r, r 2t d s , c, uj) is the coefficient depending on the lifetime of the inter- 
mediate states, T, the radial matrix element r 2 ^ S) the velocity of photon c and the incoming 



Fhki = -8y/3f (r, r 24 s, c, u) ^n ek n t > k > (t x ) [e z k z 




(4.1) 





sin Ap) + 8 sin 2 #(sin 2 p sin 2(f)] 



(4.2) 
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where 6 and ip are the Bragg and the azimuthal angles, respectively. The azimuthal-angle 
and the temperature dependences of the RXS intensities at a selected azimuth angle of the 
00 superlattice Bragg reflection (111) are shown in Fig.6. A steep decrease of the RXS 




Azimuth <P (deg.) 

FIG. 6: Temperature (Fig.a) and azimuthal dependence of RXS intensities in the ferromagnetic 
and G-type orbital ordered phase in (111) direction for the unrotated {era') (Fig.b) and the rotated 
(air 1 ) (Fig.c) channels. 

intensities at Tc is observed in Fig.6. This discontinuous decrease is associated with the lost 
of magnetic order, which diminishes a fraction of orbital interaction. Such a discontinuous 
decrease is widely observed in the RXS experiments in manganites, KCuF 3 and YV0 3 , a 
significant evidence of the strong spin-orbit coupling. The K-edge RXS peaks disappear 
completely at the OOT or MIT critical point. Therefore the RXS intensity could be a probe 
to the critical point of the MIT in orbital insulators. 

V. REMARKS AND SUMMARY 

In the preceding sections we focus on the system with FM and G-AFO structure. For 
large U/J H , the GS of the system is the A-AFM and G-AFO (C-AFO) ordered phase, as 
determined by the cluster-SCF approach. We find such an orbital ordered TMO also opens 
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an insulating gap in the single-particle energy spectra, and the energy gap exhibits simi- 
lar dependences on the 00 parameter and the temperature to those in the present paper. 
The dependence of the RXS intensity on temperature is also similar to the preceding result 
qualitatively, though the peak positions of the RXS intensities slightly shift in compari- 
son with the present structure. Therefore different ordered orbital insulators share many 
common characteristics, and are significantly different from the conventional single-band 
Mott-Hubbard insulators 1,31 \ 

Many cubic symmetric TMO with orbital degeneracy are inherently unstable, the orbital 
degeneracy is usually lifted by lowering the crystal symmetry through the JT distortions in 
LaMn03 32 and YVO3 33 , or other lattice distortion modes in V2O3 16 . As far as the JT 
phonon-mediated orbital coupling is considered, an additional orbital-orbital interaction is 
introduced 9 ); 



where g is the electron-phonon coupling constant, and K the restoring coefficient. For typ- 
ical TMO in which JT effect plays a role, the parameters g ~ 1.2 eV, and K ~ 10.0 eV. 
In this situation, the orbital correlation arises from both the Coulomb interaction and the 
JT effect, thus the orbital GS is a combination of the pure electronic spin-orbital superex- 
change interaction and the crystalline field splitting described by Eq.(2) and the JT orbital 
interaction by Eq.(14). Utilizing the cluster-SCF method, we find that the most stable GS 
magnetic structure in this case is still FM and G-AFO ordering with a slight modification 
on the orbital occupation for small U/ J H ; and for large U/ J H the GS is A-AFM and G-AFO 
ordering. Such system also exhibits most of characters of orbital insulators. 

An obvious fact is that many 00 insulators do not exhibit MIT when the long-range 00 
disappears at high temperature, such as in LaMnOs 34 -* and YVO3 33 \ except for V2O3. The 
essential reason is that these TMO are not simple orbital insulators, there also exists strong 
dynamic electron-phonon coupling, i.e., the dynamic JT effect, even if the static cooperative 
JT distortion disappears at high temperature. In these compounds, numerous dynamic JT 
phonons drag the motion of 3d electrons, and localize the 3d electrons as the incoherent 
polarons. The transport of the 3d electrons dragged by numerous dynamic phonons is of 
insulated polaronic character. Therefore these TMO are still insulators at high temperature. 
This addressed why T M ^ Too m some compounds. 




(5.1) 
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An important result from the preceding study is that the crystalline field splitting E z 
plays crucial role for the stability of the 00 GS. A pure spin-orbital interaction can not 
solely determine the magnetic and orbital GS, since the quantum fluctuation is large and 
the spin-orbital interacting system is still highly degenerate. The crystalline field splitting 
E z may suppress the quantum fluctuation and drive the system into an stable GS. Similar 
result was also obtained by Fang and Nagaosa in the orbital compound YVO3 in a recent 
paper 35 . 

In summary, in an orbital-degenerate spin-orbital interacting system, besides the orbital 
wave gap, an energy gap of the electronic excitation is opened as the orbital order develops, 
and disappears with the vanishing of the orbital order driven by the thermal fluctuations. 
The clear dependence of the insulating gap on the orbital order parameters shows that such 
TMO possesses the orbital insulator characters, and from which one could interpret the 
considerable discrepancy between the energy gap in the optical and transport experiments 
and the MIT critical temperature. The RXS intensity also exhibits the close interplay 
between spin and orbital orders. We expect more interesting properties of orbital insulator, 
such as the critical index near the transition point of the MIT, will be uncovered in the 
further studies. 
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